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Abstract. Dynamical decoupling (DD) sequences were invented to eliminate the 
direct coupling between qubit and its environment. We further investigate the 
possibility of decoupling the indirect qubit-qubit interaction induced by a common 
environment, and sucessfully find simplified solutions that preserve the bipartite 
quantum states to arbitrary order. Through analyzing the exact dynamics of the 
controlled two-qubit density matrix, we have proven that applying independent Uhrig 
Dynamical Decoupling (UDD) on each qubit will effectively eliminate both the qubit- 
environment and indirect qubit-qubit coupling to the same order as in single qubit 
case, only if orders of the two UDD sequences have different parity. More specifically, 
UDD(n) on one qubit with UDD(m) on another are able to produce min(n,m)th 
order suppression while n + m is odd. Our results can be used to reduce the pulse 
number in relevant experiments for protecting bipartite quantum states, or dynamically 
manipulate the indirect interaction within certain quantum gate and quantum bus. 
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1. Introduction 

Qubit is constantly losing its coherence due to interaction with environment. Dynamical 
decoupling sequences are proposed to eliminate the unwanted qubit-environment 
couplings with instantaneous 7r pulses, and efficiency of DD sequences is largely 
dependent on the pulse locations. Equidistant DD is a first order solution pQ, while 
more advanced sequences employ non-equidistant [21 [3] or concatenated pulse locations 
[HE]. Particularly in this article we are concerned about UDD [2], whose pulse locations 
are t5j 

5j = sin 2 (jvr/2(n + 1)). 

UDD sequence was first derived under pure dephasing models [6], able to protect the 
qubit coherence up to nth order with n pulses. 

In recent years, DD has been successfully extended to preserve multipartite 
quantum states [H El El HOl HU H21 [13] . By nesting layers of UDD (NUDD), total (n+l) m 
pulses are needed to freeze a system of m qubits to nth order without knowing any details 
of the qubit-environment coupling [H [101 HH [151 HE]- It can be seen that in NUDD pulse 
number grows polynomial with pulse order. However for a specific experimental setup, 
it is more realistic and efficient to tailor the pulse sequence accordingly, since pulse 
number is limited during a fixed time interval [T71 [18] and errors from each imperfect 
pulse will accumulate [El [201 EE]- There are already some efforts to reduce the pulse 
number. If prior knowledge is available, such as initial qubit state [7J, environmental 
coupling spectrum [U] or unbalanced decoherence rates [T5], the above universal plan 
can be greatly simplified. 

In this paper we introduce another approach to reduce the DD sequence which is 
needed for a two-qubit system linearly coupled to a common bosonic bath. This situation 
may arise when two qubits are not spatially separated enough to create independent 
environments [22], or noises felt by each qubit are correlated [Tl] 123] . In addition, 
common environment has long been exploited to generate entangling gate and serve as 
quantum bus between qubits [241 ESI ES [27]. While these gates and buses are idle we 
can use DD sequences to switch off the interactions and in the meantime protect the 
states. 

As it turns out, the simplified sequence is made very easy to implement. For 
example, we can apply UDD(n) (n-pulse UDD) on first qubit and UDD(m) on second 
qubit. n and m can be chosen at will as long as they have different parity. The whole 
sequence will preserve the bipartite state up to an order of min(n, m). For m — n + 1, 
this scheme only needs In + 1 pulses, while NUDD(n) needs (n + l) 2 . 

We organize this paper as follows. In Sec. II we discuss the free dynamics of 
a two-qubit system in a common bosonic bath. In Sec. Ill, controlled dynamics and 
conditions for high-level DD are derived. Sequences satisfying these conditions are given. 
Conclusions are put in Sec. IV. 



Dynamical Decoupling in Common Environment 



3 



2. Dynamics in Common Quantum Bath 

We consider a two-qubit spin-boson model: 

H = J2 u k b{b k + (a Zl +a Z2 )J2 h{b\ + h), (1) 

k k 

where b k is the annihilation operator for the fcth mode of the bath and X k is coupling 
strength. a zi and a Z2 are spin operators acting on first and second qubit respectively. 
This model describes a pure dephasing process due to the couplings with environment. 
With | f) and | 4) being the eigenstates of a z , we can define four basis states 
|0) = | |1) = | 41)7 1 2 ) = I 14), 1 3) = I tt) for a bipartite quantum state. 

The dynamics of Eq. (1) can be exactly solved [251128]. In the interaction picture of 
the bath operator J2k u kb\bki the unitary evolution that generated by the time-dependent 
interaction Hamiltonian can be computed using Magnus expansion. It is easy to verify 
that only the first two terms of the expansion are nonzero (see [28] for more details). 
As a result, the free evolution of the composite system can be calculated as follows 



Uf{t) = exp[(a Zl + a Z2 )Y, h( h b\) 

k ^k ^k 

, ■/ \ 2 \ 2 sin U) k t , 

+ \{a Zl + a Z2 ) 2^ X k 2 ]. (2) 

k u k 

The collective term (cr Zl +a Z2 ) 2 = 2{I+a Zl a Z2 ) generates an indirect coupling between the 
two qubits dependent on their states, which will further induce an oscillation of quantum 
correlations [251 123 121] . In order to preserve an arbitrary state, not only couplings to 
bath oscillators but also the indirect couplings must be removed. By introducing the 
spectral density function [30] defined as J{<jj) = J2k ^k^( u ~ u k), we can then write down 
the evolution of the density matrix p{t) 

I 1 e iA(t)-r(t) e iA(t)-r(t) e -4r(t) ^ 

e -iA(t)-r(t) 1 1 g -iA(*)-r(t) 



e -iA(t)-r(t) 1 1 e -iA(t)-r(t) 

„-4T(t) AA(t)-T(t) p iA(t)~T(t) -1 



with 



-r/ \ cut sin cut , 

A(t)=4/ J(u) , du, (3) 



to 2 



00 



r(t)=4/_ J(u) \ 1 / coth(^-)^. (4) 







2sin 2 (wt/2) , ,0u, 

is the inverse temperature. The exponentially decaying factor T(t) is associated with 
decoherence process, while A(t) represents collective phase evolution. Note that the 
phase factor is absent in classical noise model [HI [23]. If two qubits are subjected to 
the same classical noise, it is adequate to use simultaneous DD pulses on each qubit. 

In the case of quantum noise, still by flipping the sign of a Zl and a Z2 with tt pulses, 
T(t) will be effectively averaged to zero. However, we cannot apply the same sequence on 
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both qubits, since simultaneous n pulses will have no influence on the value of (a Zl +o~ Z2 ) 2 
and so on phase evolution. This motivates us to consider the following scenario: n pulses 
of 7r rotation along axis-x are applied to the first qubit, with the pulse locations given by 
8^,82', •••,5 n ', whereas m analogous pulses are applied to the second qubit at different 
times 8^,82", ■■■,8 m ". Totally n + m pulses are applied to the two-qubit system. We 
arrange the n + m pulse timings in increasing order and denote them by d>i, 8 2 , 8 n+m , 
with 8j < 8j + i. At each 8j, either the a zi or a Z2 operator switches its sign. At the same 
time, the operator cr Zl a Z2 changes its sign n + m times at these instants. In the next 
section we seek to find the correct n + m pulse locations which preserve both phases 
and amplitudes of the density matrix elements. 



3. Pulse Controlled Dynamics 

In this section we adopt the canonical transformation technique which was used in 
deriving the UDD sequence for single qubit [2j [6] . Also we use the notations from [6] 
by defining 

A eS = UAU\ A(t) = exp(i# eff t)Aexp(-i# eff t), (5) 
U = ex V [(a Zl +a Z2 )K],K = Y,—(bl-b k ). (6) 

k W k 

rr 

Operator A acquires time dependence under the action of the "effective" Hamiltonian. 
After the canonical transformation U, the effective Hamiltonian is diagonal 

H ee = J2u k blb k -^(a Zl +a Z2 ) 2 (7) 

k U k 

rr 

and the time-dependent flip operators <x® (t) are 

af(t) = ex P [2a Zi K(t)} exp(-4i£ ^a Zl a Z2 t)a Xv (8) 

k U k 

with 

K(t)=J2-(ble i ^ t -b k e- i ^ t ). (9) 

k U k 

Making use of these expressions, dynamics of the density matrix elements can be 
obtained by explicit calculation. For an arbitrary element p ss ' , the average evolution is 

((S\e im \S}(S'\e- [m \S')} 

e \.H(5 n+m+l -5 n+m )t ^g.^ ^ |g\ 

-lH(5 2 -h)t -\H{S t -S )t\Q'\\ 
x i\ I / / 



-eff/ 



= ((S\e 1H ^- 5 ^cr Xn e 1H ^- 


-5i)t ^lH(6 n+m - 


-£ n _|_ m _l)£ 

U X 


Jn-\-m 


(S' 1 (a x .)e- iH{5n+m+1 ' Sn+m)t a x . e -' m& ^ 


■m ^n+m — 1 )^ 


...e~ 






\-m-lt)o~ Xj ^ 


(* 

-m x 


(S'\U^ iHe \f(t))a x ^ 


■m v 1 / Jn- 


fm-i \^n+m- 


-it). 




■ ■■o- x (8 n - 


\-m-l't) 'x jn _ t 


(* 

-m x 
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e(--)^)e- 1 ^^ (si2 -- 2) V)^lKf(0)^„ + JW) 

&X4 



, ,„ {5 n + m -it)...<J Xj2 (S 2 t)a Xji (6!t)U\S )>. 

Si, S2 are eigenvalues of cr^ + <r Z2 for IS*) and IS''). Determined by which spin is flipped 
at 5j, ji, J2, ...jn+m take values from {1,2}. By (o~ Xi ) we mean that an ending pulse a Xi 
is added at the end of the sequence to ensure the final state is still |S)(S | if n + m is 
odd (If n and m are both odd, two ending pulses are needed). However according to 
our calculation, whether or not there is an (or two) ending pulse will not modify the 
equations and results we derive next. The outermost brackets denote ensemble average 
with respect to the thermal bath. 

The sum of all terms in the form of — 4i Y,k{^\l UJ k)o'z 1 o' Z2 from Eq. (10) leads to a 
phase shift 

exp[-4i]T - S2 )(5 X - 5 2 + ... + {-l) n+m - l 5 n+m + ^-L )t], 

k Uk 1 

(11) 

and calculation of the rest part of Eq. (10) will produce more phases. Taking poi(0 as 
an example, we are left with 

^ e 2K e -2K(5 1 t) e (±)2K(5 2 t)^J±)2K(S rn+n t) e (-l)^2K(t) 
e (±)2K(S m+n t) ^ e (±)2K(S 2 t) e -2K(S 1 t)y 

(12) 

Coefficients before each K(Sit) are determined by the relative locations between lst- 
qubit and 2nd-qubit pulses in the whole sequence. Exponentials in Eq. (12) can be 
combined one by one using Baker-Campbell-Hausdorff formula e A e B = e A + B e l A < B }/ 2 ^ 
which is valid here because 

[K(t),K(t)} = 2^1/^)^(1-1) 

k 

is a c-number. After the combination all K(5it) add up to one exponential 
exp[-2^(Vt) + 2K{5 2 >d)... + {-l) m 2K{5 m/ ,t) + (-l) m+l 2K(t) + (-l) m 2K(<5 m »t)... + 
2K(5 2 "t) — 2K(5 1 "t)] while extra phases introduced by the e^ A ' 3 ^ 2 term in Hausdorff 
formula can be calculated by observing that combining any four symmetric exponentials 
e ±2K(8 p t) an( j e ±2K(8 q t) on both sides of e (- 1 ) m+l2A: (*) w ill not create extra phases except 
for one case: that is when combine e^ -1 ^ 2K ^ i ' t "> with the exponentials which act on the 
2nd qubit after it 

e (-iy2K(S.,t) e (-iy2K(S / ,t) e (-iy2K(S / n) e{ - 1 y { -2)K(S. l t) 



(-i)*+n6i£ fc J| sin Wfc (5 ,-5 , /)t ^y^^,,^ , .„ 
; fc J e i ,(i < j ). 



(13) 

Merging these exponentials four by four, we arrive at 

\ 2 n 

ex P {i£%8E(-l) i+m+1 sm^' _ l) t 

k U k j=i 
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+ 16E £(-l)^'sinu; fc (^-<V)*]} 



7 — 1 •' '" 
l — 1 I <J 



(exp(2AT) exp[-2K(5 1 »t) + 2K(6 2 >d)... 

+ (-l) m 2K(5 m „t) + (-l) m+1 2K(t) + {-lT2K{5 m>l t) 

... + 2K(5 2 »t) -2K(8 1 »t)\). (14) 

The ensemble average in Eq. (14) has already been calculated in [6]. With the existing 
results from [6] and combining Eq. (11) and Eq. (14), the final form of poi(t) can be 
organized in a rather compact way 

p 01 (t) = exp(iA m (i) - T m (t)), (15) 



with 



roo J((jj\ 

A m (t) =4/ — —{x m + z m + c)(ut)doj, 



UT 

oo 



\: ni (t ) 2/ J(u) \ Vm ^ 2 coth(pu/2)du, 



U z 



x m (ut) = (-l) m sin(u;i) + 2 £(-l) i+1 sin(w%), 



</ m (u;t) =l + (-rr+V"' + 2£(-l)^ 



7 lwt(5 // 



^(w*) =2X:(-l) l+m+1 sin^-l) + 4£ ^(-l^sinwt^-^), 

i=\ *=i i'</' 

n+m / i \n+m 

c(wt) = -2o;t[£(-l)'- 1 5 r + ^ ]. (16) 

Unlike single qubit DD, here the phase factor A m (t) will induce collective dynamics and 
should be minimized. y m {uii) is the so-called filter function of an m-pulse sequence [21 Hj 
and T m (t) is responsible for decoherence. 

The most simple way to realize dynamical decoupling is to use independent UDD 
sequences to suppress the decoherence, and hope that phase evolutions will be minimized 
automatically at the same time. UDD(m) requires the first m derivatives of the filter 
function to be zero at uj = 0. These constraints are imposed to engineer y m (ut) in low 
frequency region, achieving an error order of 0(t m ). Similarly we can also define the 
filter function for A m (t) as 

f m (ut) = X m (ut) + Z m (ut) + c(u)t). (17) 

The density matrix p(t) that controlled by an arbitrary pulse sequence reads 

/ i e iA m (t)-r m (i) e iA„(t)-r„(t) e -r„(t)-r m (t)- 7 (t) \ 

e -iA m (t)-r m (t) ]_ e -r«(t)-r m (*)-7(t) e -iA n (t)-r n (t) 

e -iA n (i)-r„(t) e -r„(t)-r m (t)-7(t) ^ e -iA m (t)-r m (t) 

i e -r„(t)-r m (t)- 7 (t) e iA n (t)-r n (t) giA TO (<)-r TO (t) -j_ 
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Figure 1. t = 1. Red line: UDD(2) on the first qubit, UDD(3) on the second 
qubit. Solid line for / 2 , and dashed line for / 3 . Other four colored lines follow similar 
definitions. 

with 

7 (t) = 4 ^.^(y n (u;t)y* m (ut)) coth(^/2)dw. (18) 

JO UJ Z 

Definitions of A ra (t) and T n (t) are the same as Eq. (16), but with pulse locations 
swapped. For example, 

m 

z n (cot) = 2]T(-iy +n + 1 smu;t(<5 i » - 1) 

3=1 
m n 

+ 4 E E(-!) i+isin ^%'-^)- (19) 

P12 and P21 are driven out of the decoherence-free subspace. In spite of this, it is clear 
that the exponential decay of all density matrix elements are filtered by y m (ut) and 
y n (uit), which is at the order of min(n,m) if UDD(n) and UDD(m) are applied on each 
qubit. 

3.1. UDD Sequences with Different Parity 

By requiring the phase shift from Eq. (11) vanish, we get the equality 

/ 1 \n+m 

6,-62 + ...{-l) n+m - l 6 n+m + ^-L = 0. (20) 

This is exactly the equation to derive UDD(l). In other words, our (n + m)-pulse 
sequence has to be a first order sequence at least. Now we make the first observation: 
Eq. (20) holds for any combination of UDD(n) and UDD(m) with n + m odd. Besides, 
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parity plays an important role in the following relations 

m n 

E E(-l) t+ ^m^(V-^) 

n m 

= EE (-l) 4+J sinwt(5 (n+1 „ i) / -5 (m+1 _ jr ), 

i=li'>j" 



and 



= E E (-l) i+i+m+n sinu,*(<V -8 f ), (21) 

* <3 
m 

J2(-iy +n+1 smtut(5 f - 1) 

3=1 

m 

= -E(-l) J+rt+1 sm^(5 (m+1 _, r ), 
i=i 

m 

= -£(-l)J +n+m sina;t(fy/). (22) 

As a result, if n and m have different parity, we have 

x m {ujt) + z m (ut) + x n (ut) + z„(u;£) = 0. (23) 
Moreover, it is easy to verify 

x m {ut) + z m (ut) - x n (ut) - z n (wt) = %(y n (t)y* m (t)), 

(24) 



since 



£ (-1)*+* Bmutfy - 5, 



7 — 1 ^ -ff 



- E E(- 1 ) i+ ' sinw %'-^') 

n m 

= E E(-ir +J smc^-V) 

i— 1 ■' ^ " 
m n 

+ E E(- 1 )' +isin ^0V-<V')> 
1—1 " ^ ■' 

n m 

= EE^ 1 )^^^-^)- (25) 
i=i j=i 

For n + m is odd, solving Eq. (23) and Eq. (24) yields 

f m {ut) = -f n (ut) = <&(y n (t)y* m (t))/2. (26) 

Again the phase factors A m (i) and A n (t) are bounded by the filter functions y m (t) 
and y n (t). Thus we have completed the proof that both phase and coherence can be 
preserved to the same order. In Fig. [TJ we give three examples of independent UDD 
sequences with different parity, where f m {uit) and f n (ut) are numerically calculated. It 
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Figure 2. The same definition as Fig. [TJ Note that the solid lines completely overlap 
the dashed counterparts as indicated by Eq. (27). 

can be seen that the phases are suppressed order by order through increasing the pulse 
number of each UDD sequence. 

3.2. UDD Sequences with Same Parity 

We only consider two UDD sequences of even order, and two odd-order sequences share 
essentially the same property because their middle pulses do not flip the collective 
operator cr Zl o~ Z2 . If n and m are even, from Eq. (21) and Eq. (22) we know 

x n {ut) + z n (ut) = x m (ut) + z m (ut). (27) 

A corollary is drawn: y n (t)y* m (t) is a real number if n and m are even. Moreover, 
x m {ujt) + z m (ut) are not bounded by filter functions anymore. As shown in Fig. |2j 
phase evolutions are eliminated at a fixed level. Increasing pulse number cannot improve 
the performance. Suppression of f n (ut) and f m (ut) before ut < 2 is due to the fact: 
sincjt ~ cut if ut is small. As a consequence, f n (ut) and f m {ut) are self-corrected to 
first order, see Eq. (2). While a zl a Z2 ut is effectively averaged at low frequencies, the 
nonlinear part of cr Zl a Z2 sin cut is uncontrollable. 

In the present model the indirect coupling via a thermal bath is commonly weak. 
If there exists highly nonlinear indirect coupling between the qubits, we expect more 
distinctions will be observed by using different UDD sequences. 

3.3. Entanglement Dynamics 

In order to illustrate the difference caused by the parity, especially the distinct oscillation 
patterns of quantum correlation under different decoupling schemes, we numerically 
calculate the entanglement dynamics. The initial bipartite state is chosen to be 
^(|0) + a/3|1) + \/3|2) + |3)), which is partially entangled. In Fig. [3] we have used 
concurrence to measure the two-qubit entanglement [31 J . 
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Figure 3. Evolutions of concurrence in three different environments. O(-) is Heaviside 
step function. In each of the three diagrams, red line stands for the same parity, and 
the other colored line for different parity. 

We consider the dynamical evolution of concurrence versus time under three kinds of 
environmental spectrum. We can see that for J (to) = to 0(1 — to), the final concurrence 
at t = 1 is perfectly kept at the initial level regardless of the parity. This result is 
consistent with our earlier observation that the low frequency part of f(tot) is self- 
corrected. However for the other two spectrums, the UDD performances are obviously 
dependent on the parity difference. For J (to) = toO(5 — to), the combination of UDD (6) 
and UDD (7) still preserves the initial entanglement at the end of the decoupling cycle, 
while UDD(6) and UDD(8) fails to do so. Since UDD(6) and UDD(8) are not able 
to effectively suppress the indirect coupling between the two qubits, the quantum 
correlation of their final state is much larger than the initial one. Similar increase can 
be observed in soft-cutoff spectrum (J (to) = j^p)- Note that in this type of spectrum 
concurrences cannot be well preserved in both combinations owing to severe decoherence 

©EE]. 

As shown in Fig. [3j concurrences undergo violent oscillations within the pulse 
intervals. In contrast with UDD(6) and UDD(8), the concurrence that is controlled 
by the combination of UDD(6) and UDD(7) drops more rapidly in the beginning. For 
Ji 00 ) = i+b 7 ' ^ ne concurrence even reach zero at one point, indicating the entanglement 
is completely lost. However at the final stage of the decoupling period, the sequences 
with different parity have a more smooth and steady concurrence, which makes them 
more feasible since the loose timing constraint for retrieving the final state, compared 
with the combination of UDD (6) and UDD (8). 

4. Conclusion 

In this paper we give a detailed analysis of available DD sequences which can be used 
to eliminate both qubit-environment and indirect qubit-qubit coupling in a common 



Dynamical Decoupling in Common Environment 



11 



bath. Exact dynamics under arbitrary pulse sequences are derived. As we find out, it 
is possible to apply UDD sequences independently on each of the two qubit, and at the 
same time preserve the bipartite state to higher order. This result greatly simplifies 
the scheme of universal NUDD when dealing with correlated environments. We has 
proven that by applying UDD(n) and UDD(m) with n + m odd, the evolutions of all 
density matrix elements are bounded by the single-qubit filter functions of UDD(n) and 
UDD(m). 

In conclusion, we suggest using UDD sequences with different parity to protect 
a dephasing bipartite system, in case that the environments for individual qubit are 
not completely independent. Besides, our results may find applications in quantum 
information processing due to its superior ability to dynamically switch off the 
interaction induced by a common quantum bus. 
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